Spin-orbital effects in magnetized quantum wires and spin chains 
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We present analysis of the interacting quantum wire problem in the presence of magnetic field 
and spin-orbital interaction. We show that an interesting interplay of Zeeman and spin-orbit terms, 
facilitated by the electron-electron interaction, results in the spin-density wave (SDW) state when 
the magnetic field and spin-orbit axes are orthogonal. We show that this instability is enhanced in 
a closely related problem of Heisenberg spin chain with asymmetric uniform Dzyaloshinskii-Moriya 
(DM) interaction. Magnetic field in the direction perpendicular to the DM anisotropy axis results in 
staggered long-range magnetic order along the orthogonal to the applied field direction. We explore 
consequences of the uniform DM interaction for the electron spin resonance (ESR) measurements, 
and point out that they provide way to probe right- and left-moving excitations of the spin chain 
separately. 

PACS numbers: 71.70.Ej,73.63.Nm,75.10.Pq 
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I. INTRODUCTION 

Over the last several years there has been a remark- 
able growth in research activity in the field of spintron- 
ics with the ultimate goal to fabricate novel spin-filter 
devices which can control and manipulate the electron 
spinsi. The proposals for such a spin-filter device con- 
tains two achievable attributes, a ballistic quantum wire 
and the presence of a tunable Rashba spin-orbit coupling 
responsible for controlling the electron spin. Ballistic 
quantum wires are created in a 2DEG by cleaved edge 
over growth, whereas the Rashba effect arises due to the 
asymmetry associated with the confinement potential 2 . 
The asymmetry and hence the Rashba coupling strength 
can be controlled by applying the gate voltage. Al- 
though the role of spin-orbital and magnetic (Zeeman) 
fields in the electric and spin transport is well under- 
stood for a non-interacting quantum wire^££i£, the 
case of interacting electrons remains the subject of ac- 
tive TeseardM^^M^^e. 

It should be noted that finite spin-orbit coupling is 
very natural, and, strictly speaking, unavoidable, in semi- 
conducting quantum wires due to pronounced structural 
asymmetry inherent in the fabrication process. Also, in 
addition to the noted asymmetry of confining potentials 
(which include quantum-well potential that confines elec- 
trons to a 2D layer as well as transverse [in-plane] po- 
tential that forms the one-dimensional channel 3 -), spin- 
orbit interaction is inherent to semiconductors of either 
zinc-blende or wurtzite lattice structures lacking inver- 
sion symmetry^. 

Another very interesting system that motivates our in- 
vestigation is provided by one-dimensional electron sur- 
face states on vicinal surface of goldi 8 - as well as by 
electron states of self-assembled gold chains on stepped 
Si(lll) surface of silicon-!-?. In both of the systems 
one-dimensional ballistic channels appear due to atomic 
reconstruction of surface layer of atoms, see also 20 ' 21 . 
The resultant surface electronic states lie within the 
bandgap of bulk states, and thus, to high accuracy, 



are decoupled from electrons in the bulk. Spin-orbit 
interaction is unexpectedly strong in these systems, 
with the spin-orbit energy splitting of the order of 
100 meV. In fact, spin-split subbands of Rashba type 
have been observed in angular resolved photoemission 
spectroscopy (ARPES) in both two-dimensional 22 and 
one-dimensional setting a 18 ' 19 . The very fact that the 
two (horizontally) spin-split parabolas are observed in 
ARPES speaks for high quality and periodicity of the 
obtained surface channels. 

As we show below, the most interesting situation in- 
volves electrons subjected to both spin-orbital and mag- 
netic fields. While it is perhaps impractical to think of 
ARPES measurements in the presence of magnetic field, 
it is quite possible to imagine experiments on magnetic 
metal surface a 23 ' 24 . It is then natural to investigate com- 
bined effect of non-commuting spin-orbit and Zeeman in- 
teractions, together with electron-electron interaction, on 
the one-dimensional system of electrons. 

Electrons in a quantum wire (or, in a one-dimensional 
surface channel) are a good realization for a Tomonaga- 
Luttinger liquid and serves as an ideal system for the 
study of the interplay of magnetic field and Rashba 
spin-orbit effect on the interacting quantum wire. The 
magnetic field breaks the time-reversal symmetry of the 
Hamiltonian and splits the band of free electrons into 
two, corresponding to up-spin and down-spin electrons, 
reducing spin-rotational symmetry of the system from 
SU(2) to U(l). Subsequent inclusion of the Rashba term, 
Hr oc a X p- z, see Eq.©, breaks this U(l) symmetry 
(observe that [<r z ,Hn] ^ 0). In addition, the spin-orbit 
interaction (SOI) Hr breaks spatial inversion symme- 
try V : x — > — x. A consequence of the fully broken 
SU(2) symmetry is the generation of new scattering pro- 
cesses which are no longer spin-conserving. These are 
the Cooper scattering processes in which a pair of elec- 
trons in the lower band scatter to the upper band and 
vice versa££. A relevant Cooper term creates a gap in 
the energy spectrum leading to a long range spin-density 
wave order. We find that the gap strength is proportional 
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to the backscattering (2k f component) of the electron- 
electron interaction potential and the ratio of the Rashba 
to the Zeeman energy. For a large enough gap the order- 
ing in the spin-density wave can crucially suppress the 
backscattering process of electrons from an isolated im- 
purity. Brief description of our main results was previ- 
ously given in Ref . [25l 

We will also analyze an alternate system, a Mott- 
Hubbard Heisenberg spin-1/2 chain, in the presence of 
magnetic field and Dzyloshinskii-Moriya (DM) interac- 
tion. A magnetic field applied at an angle perpendic- 
ular to the DM anisotropy axis breaks the continuous 
U(l) symmetry and consequently a true long range or- 
der can develop in the spin-chain. The case of a stag- 
gered DM term has been studied experimentally — and 
theoretically 27 , and has been shown to open up a gap 
in the spectrum with the gap scaling as B 2 / 3 with the 
magnetic field B. The case of a uniform DM term is 
also experimentally relevant, see for example Refl28l. but 
has not been discussed much theoretically. We will show 
that the case of a uniform DM term and perpendicular 
magnetic field can be described analogously to the quan- 
tum wire in the presence of spin-orbit interaction and 
magnetic field. 

The outline for the paper is as follows: In Sec|TT] we 
will review non-interacting electrons in the presence of 
magnetic field and Rashba spin-orbit term. We then 
consider interaction effects by using standard bosoniza- 
tion approach. In Sec lIIIBl we will perform renormaliza- 
tion group analysis to determine relevant and irrelevant 
terms. In SeclIVI we use perturbative approach to gen- 
erate relevant terms. The role of relevant terms on the 
transport property of electrons in the presence of a sin- 
gle impurity is analyzed in section [Vj In Scc lVIi we con- 
sider a spin-1/2 Heisenberg antiferromagnetic chain in 
the presence of magnetic field and uniform DM term and 
compare this system with the quantum wire. Compari- 
son with the case of a staggered DM term is considered 
in Sec I VI Dl In Section IVIII we discuss the role ESR mea- 
surements can play in unraveling the DM term in the spin 
chain. Technical details of our calculations are described 
in Appendices. 



(fx = x, y, z) are the Pauli matrices and the potential 
V(x) = muj 2 x 2 j2 typically confines the particle in the 
x— direction. When the confining potential is strong 
enough so that the width of the wire \J hj(2muS) is much 
smaller than the electron Fermi- wavelength, only the first 
sub-band is occupied and the Hamiltonian {T]) acquires a 
one-dimensional form 



Hn 



h 2 k 2 

2m 



(a x k) 



a 

■ B. 



(3) 



Here and below k is electron's momentum along the axis 
of the wire, which we will denote as x-axis in the following 
for notational convenience. We also set ft = 1. 

It is easy to see that in the absence of magnetic field 
SOI in © can be easily gauged away via the spin- 
dependent shift of the momentum, H (B = 0) cx (k + 
ma.RO x ) 2 . Corrections arising from the omitted term 
a RCyPx © produce small spin-dependent variations of 
the velocities of right- and left-moving particles^. These, 
however, are not important for our purposes for as long as 
ctRkp <C Ep, which is the limit (along with Az <C Ep) 
considered in this work. With interactions included, elec- 
trons form Luttinger liquid with somewhat modified crit- 
ical exponents, in comparison with the standard case of 
no SOli&ii. 

Most interesting situation arises when both SOI and 
Zeeman terms are present simultaneously and do not 
commute with each other as happens when spin-orbital 
axis (a x in ([3])) is different from the magnetic field direc- 
tion. In what follows we choose magnetic field to point 
along the z— direction, B = Bz. The energy eigenvalues 
of the Hamiltonian j3|) is found as^ 



e± = ^±J(aRk) 2 + (^) 2 , 
2m V 2 



(4) 



and the momentum dependent eigen-spinors are, for 
e+(k): 



\x + (k)) 



sm 



7(fc) 
2 

7(fc) 



(5) 



II. ID ELECTRONS IN THE PRESENCE OF 
MAGNETIC FIELD AND RASHBA SPIN ORBIT 
TERM 



and for e-(k): 



|X-(*0) 



cos 



— sm 



7(fc) 
2 

7(fc) 



(6) 



The Hamiltonian for an electron subject to the Rashba 
spin-orbit term and in the presence of magnetic field is 
given by 



Ho = h(Px 2 * Py) +V(x)- 9 ^.B + Hr, (1) 

(2) 



ii aR ( \ 

H-R = -^-(<J x Py - <TyP x ) 



where aR is the Rashba spin-obrit coupling, g is the 
effective Bohr magneton, B is the magnetic field, 



where A z — gpsB and rotation angle j(k) is introduced 



7(fc) = arctan 



2a R k 
A z • 



(7) 



Notice that particles with momentum ±fc experiences an 
effective magnetic field B e g = ^2aRk/(gpB)x + Bz, thus 
spin directions in each band varies with the momentum. 
Going from left to the right side of the band, the spins 
"rotate" along the counter-clockwise direction. At k = 0, 
the effective magnetic field is just the applied field, thus 
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the separation between the two bands is minimum and 
the spins align according to the applied field, i.e., along 
the =Fz— direction. For states with the same energy the 
spin states of + and — bands are no longer orthogonal 
if there is a finite magnetic field and Rashba spin-orbit 
coupling. In particular the right and left Fermi levels 
satisfy the following property 



k R/L + kl ,L 



±2k f 



and 



5k F = \k R/L -kl /L \ 



■^/(2k F a R ) 2 + (A z 



The magnitude of velocities at the Fermi level are 



v F =F 



(A, 



2k F y / A(a R k I 



(8) 



(9) 



(10) 



where v F = k F /m. The spin overlap between the upper 
(+) and lower (— ) band is non-zero, 



(x(V )\x(kJ ))=sm ± . (11) 

As will be discussed in the next section, the non- 
orthogonality of the spin-states acquires important con- 
sequences when one turns on the electron-electron inter- 
actions. 




^> a (x) is conventionally defined in terms of the annihi- 
lation operator a a (k) of a free particle in the state e tkx 
and with spin a =t,|: 4^0) = / ^e lkx a a {k). Alter- 
natively, annihilation operators a±{k) of particles which 
are the eigen-states of the Hamiltonian {T]) with eigen- 
energies e±(k), can be used to represent the field operator 
as follows: 



u=± 



dk 
2^ 



e ikx { Xv {k)\cj)a v (k). (13) 



The low energy physics of the interacting wire is de- 
scribed by linearizing the spectrum near the Fermi- 
points, ±k±. The a =\,\ field operators are now de- 
scribed, in coordinate space, in terms of the chiral right 
(Rv) and left (L v ) movers of v — qp subbands as follows 



{ Xv {k v )\a)e^ x R v + { Xv {-k v )\o)er^ L 



Following^-, we decompose the interaction part of the 
Hamiltonian (|12[) into intra-subband ifintra and inter- 
subband iJ; n ter scattering processes. The intra-subband 
process, H ln t r& , describes the interaction between elec- 
trons in the same subband and involves the standard for- 
ward and backscattering processes, see Appendix [Al The 
second scattering mechanism, -ffintcr, involves scattering 
between electrons in different subbands and can be conve- 
niently divided into forward, backward and Cooper scat- 
tering processes. Below we will discuss the inter-subband 
scattering processes in more detail. The forward scatter- 
ing process involves interaction between q compo- 
nents of the densities in the two subbands 

Hl tel = Jt/(0) f dx (RlR» + LlLv) 



x{Rl v R- v + I7_ V L- V ). 



(15) 



FIG. 1: Occupied subbands e± of Eq. [4] Arrows illustrate 
spin polarization in different subbands. Dashed (dotted) lines 
indicate exchange (direct) Cooper scattering processes. 



III. INTERSUBBAND INTERACTION 
EFFECTS 

The interaction part in terms of the particle-field op- 
erators ^ a (x) and tyJrix) ( a and a' are the spin indices) 
is given by 

Hi*.t = \ J dxdx'Uix-x'^Kx^lix'^^ix'^^x), 

(12) 

where the summation on pairs of identical spin indices 
is assumed and U(x — x') is the screened (by surround- 
ing gates) interaction between the electrons. The field 



The inter-subband backscattering process is classified 
into direct and exchange scattering. Direct backscatter- 
ing process involves q ~ 2k± components of the densities 
in the two subbands: a left (right) moving electron in the 
subband v (—v) changes its direction to become a right 
(left) moving one while remaining in the same band, 



d-B 



inter 



= cos[7(fc + )] cos[7(fc_)]C/(fc + + fc_) 



[ dx e i2 ^- k -^ x (RtL u )(Lt v R. u ). (16) 



This contribution involves an oscillatory factor 
exp[i2(fc„ — k- u )x] in the integral due to the non- 
conservation of momentum during the scattering. 

The other backscattering process is via a momentum 
conserving exchange mechanism, where electrons again 
scatter by large momentum transfer and in the process 
exchange their bands. The scattering channel conserves 
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momentum and reads 



MK) -i{k- v ). 



v=± 



x / dx^RiR- u )(RlM + {LlL- v ){l)_ v L v j) 
+2U{k lJ + fc_„)sin 2 



2r7(M+7(fc-r/) 1 



x / dx 



(RiL_ v )(LlM}. 



(17) 



Note the appearance of (squared) wave function overlap 
factors, ex (x(&j,)|x(±fc-i;))> which signify the exchange 
nature of the scattering. Note also that these factors 
are non-zero due to a finite Rasbha coupling an, which 
allows electrons to scatter without conserving their spins. 

The Cooper scattering process, which is central to our 
story, involves scattering of a pair of opposite movers 
(right and left) in the subband v into a similar pair in 
the other, — v, subband. Each pair has zero total mo- 
mentum which remains conserved in this scattering. Be- 
ing pair-tunneling like, Cooper scattering requires non- 
conservation of spin. It represents, for example, a scat- 
tering of a pair of two (almost) down-spin electrons into a 
pair of two (almost) up-spin ones. The Cooper scattering 
reads 



{u(k- 



c CI = jd X {u(k -/,,.sur- (/ '' "'I 

U(k- + fc+)sin 2r 
?t rt 



r 7(fc_)+ 7 (fc+) 



1} 



x(R 1 _L 1 _R+L + +h.c. 



(18) 



The first term (direct Cooper scattering) is due to elec- 
trons in the band R v and L v jumping into and 
L- v , respectively. The coefficient for this term is oc 
U(Skp) sin 2 [(7_ — j + )/2], where Sk F = \k v — k^ v \ is the 
momentum transfer for an electron and sin 2 [(7_ — 7 + )/2] 
is the squared overlap integral. (For brevity, we de- 
note 7(fe±) = 7± here and in the following.) The sec- 
ond process (exchange) with electron scattering from R v 
and L„ to L^ v and R~ v , respectively, involves a coef- 
ficient, cx U(k u + fc_^)sin 2 [(7_ + 7 + )/2], with a larger, 
sin 2 [(7_ + 7 + )/2], overlap integral. The bigger overlap 
for this second (R v <-* L- u ) process is also rather clear 
from pictorial representation of spin orientation in differ- 
ent subbands, as shown in Fig. [1] For the case of short- 
ranged (screened) interaction potential a simple estimate, 
using © and ©, 



U(k- 



k + ) sm 



2 [7- ~7+ 1 



U(k- + k+) sin' 



2r7-+7+ - 



U(Sk f 
U(2k F ) 



(M « > <»> 



shows that the second, exchange Cooper process, domi- 
nates. This defines the regime to be considered in this 
work. 

Finally we take into account two classes of momentum 
non-conserving scattering processes where one of them 




FIG. 2: Asymmetric back-scattering processes. 



exhibits mixed features of Cooper and back-scattering 
(called as asymmetric back-scattering process) while the 
other one has features reminiscent of Cooper and forward 
scattering processes ( asymmetric forward-scattering pro- 
cess). A typical asymmetric back (forward) scattering 
event involves total momentum change ±6k F , see FigfS] 
For example, right and left moving fermions in the same 
sub-band scatter into left (right) and right (left) moving 
fermions respectively, with one of the fermions now in 
a different sub-band. Alternatively, the oppositely mov- 
ing fermions may be in different sub-bands to begin with 
but end up in the same sub-band with opposite (same) 
momentums. The scattering process acquires a slowly 
(compared to direct back scattering (fl7j|) ) oscillating fac- 
tor exp[i(fc^ — k- u )x] owing to the non-conservation of 
momentum and is given by 



H, 



asymm 



sin( 



7+ + 7- 



^2 sga(v)U(2k F ) cos(7„) 



dx[e l ^- k -^ x {RlLlR v L^ - R}_ v L\R v L v ) + h.c.](20) 



The above expression reflects contributions from only the 
asymmetric back-scattering processes. The asymmetric 
forward-scattering processes involve identical fermion op- 
erators as in Eq. (f2"D)) . However the ratio of amplitudes for 
asymmetric forward to asymmetric back-scattering pro- 
cess is small, 



U(5k 



F sin 



U(Sk F ) y/Al + A(a R k F y 



U{2k F )sm 



■TV+7-j. 



]cos7„ U(2k F ) 



2E F 



< 1, 



which allows us to neglect contributions from asymmetric 
forward-scattering processes altogether. 

The electron density in the quantum wire is assumed 
to be incommensurate with the lattice spacing, hence the 
Umklapp scattering process is not considered. To sum- 
marize, the interaction part of the Hamiltonian has been 
decomposed in terms of three broadly defined scattering 
processes, intra subband, inter subband and asymmetric 
scattering process, 



H — H u 



(21) 
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A. Bosonization 

Bosonization is performed by expressing the fermionic 
operators in the Hamiltonian via the chiral bosonic <Pr/l u 
field s 30 ! 31 ' 32 . The fermionic fields in terms of the chiral 
bosonic field are as follows 

R± = V± p'^<^± L±= V± P -'^*-± ,(22) 
\/2™o \/2™o 

where ao ~ kp 1 is the short distance cutoff and r]± are the 
Klein factors which are introduced to ensure the correct 
anti-commutation relations for the fermionic operators 
from different (±) subbands. The bosonic operators obey 
the following commutation relations: 



f>Rv,<t>Lu>] = j<W; where i/,i/ = ± (23) 



[(t>R/Lu{x),<t>R/L V > (y)} = ±-(Wsign(a; - y), (24) 

the first of which, (|23[) . ensures anticommutation between 
right and left movers from the same subband, while the 
second is needed for the anticommutation between like 
species (i.e. right with right, left with left). Klein factors 
anticommute 



{t)»,V} = 2< W j vt = Vv 



(25) 



In the following we choose the gauge where r\-\-r\- = i. 
The chiral 4>n/Lv are expressed in terms of <j> v and its 
dual as follows 



9Rv 



"Lis = " ; 



2 ' 



(26) 



The bosonized form of the Hamiltonian is obtained by 
making use of equations (f2"2"|) through (f2"6")l , as well as the 
following results for (chiral) densities 



R v Ru 



dx<t>Rv dx(<Pu - V ) 



4tt 



Li Lij = ^ = ^ + M (2?) 



4-7T 

The (bosonized) Hamiltonian in terms of <p v and Q v is 
the sum of intra-subband, inter-subband scattering and 
asymmetric scattering processes, 



H ^intra ~\~ "inter ~t~ -"asymm. 



(28) 



where the intra-subband part has the usual form (see 
Appendix |A"|) 



flintra = \ ^ [ dx{v F (d x 9 u ) 2 

+ ( j- ^(0)-cos 2 [ 7 (fc,)]t/(2fc,) wa 2l 

+ («fH ){o x (f>u) f-(29) 

7T J 



The inter-subband part of the Hamiltonian, Hi nter , is 
given by 



tt ttF I rrd-B , rr( 

Winter — -"int.dr I "int.Ar T "i 



:-B 



inter — ■« inter "T" n inter "1" -"inter "T" ^inter, 



(30) 



where 



m 



U(0) 



dx d x (j) + d x (j)- 



d _ B {/(2fc F )cos 2 [ 7F ] 
inter 2(™ ) 2 

x / c?xcos[2v / 7r(</)- — 0+) + 2(fc + — fc_)a:], 



HI 



cx-B 



E/(2fc F ) . 



2tt 



■sm [7fJ 



x dx {d x (j) + d x 4>- - d x 6 + d x 6-), 



H 



G l/(2fc F )sin 2 [ 7F ] 
inter 2(7ra ) 2 

y dxcos[2V7F(^- - 



(31) 



The asymmetric part has the following bosonized form, 



H. 



V2U(2k F ) sin(27 F ) 



(27r) 3 / 2 a 

x sin[v / 47r(</)L- - <Pl+) - Sk F x] + d x (<p L 
x sin[V47r(</)_ R _ - <j) R+ ) - Sk F x}}. (32) 



dx{d x (4> R - - 4> R+ ) 



In deriving Eqs. (|3"Tj) and (j3"2")l we took limits A 2 <§; i? F = 
v F k F and a R kp <C Ep which allowed us to neglect veloc- 
ity differences (|10p in the two subbands and approximate 
U(2k ± ),U{ * k »- k -" ) « t/(2fc F ) and 7 (fc±) » 7 (fc F ) = 
7 F . An important exception to this replacement is pro- 
vided by -fTitter and #asymm in (JSJ) and (021), respectively, 
where momentum mismatch factors 2(fc_ — fc + )x = 2Skpx 
and (fc_ — fc + )x = Skpx must be preserved. It is worth 
noting here that the approximations assumed do not re- 
strict the ratio 2a R k F /(gfisB) = E s _ a /A Zl which can 
still take on any value. 

A more standard representation of the Hamiltonian 
is in terms of the symmetric <p p , 9 p (charge) and anti- 
symmetric (p^, 6 a (spin) modes. These combinations are 
defined as follows 



V2 ' 



V2 



V2 ' 

The Hamiltonian now reads 



V2 



H — Hp + H a 



(33) 



(34) 



The charge part of the Hamiltonian is harmonic 



H p = - I d.r 



u p Kp(d x p 



K 



^(d x <p p ) 2 . (35) 
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The spin part is the sum of quadratic and non-linear 
terms, H a = H° + H° + Hf + , where 



Hi 
H° 
H* 



\ 1 



C/(2fc F )sin 2 [ 7F ] 

2(™ ) 2 
t/(2fc F )cos 2 [ 7F ] 

2(™ ) 2 

^asymm — ij]-\-7]- 

dx{d x ((j) 



(36) 
(37) 

dxcos[V8nip a + 2Skpx] (38) 
(7(2fc F )sin(2 7F ) 



J dx coslV&rrOcr] 



(27T)3/2a 

9 a ) sin[V2^((/ )(T + 6> ff ) - Sk F x] 

+d x {^ a + 9a) sin[\/2^ CT - 6 a ) - 5k F x}}. (39) 

For completeness, it is worth noting that the leading cor- 
rection to these equations is represented by the inter- 
mode term 

H'p-o = (M+ ~ UJ) [dx{d x <p p d x <p a + d x 6 p d x d„), (40) 

which couples spin and charge sectors. Its small ampli- 
tude (u+ — U-) oc A z /Ef, see ([10]). justifies its neglect 
in the following. 

Competing nature of interacting problem is clear from 
the presence of two non- linear terms, (|37|) and (|38|) . in- 
volving non-commuting (dual) boson fields 9 a and ip a , in 
the Hamiltonian. Similar situation happens in models of 
organic conductors, where spin-nonconserving spin-orbit 
and dipole-dipole interactions play an important role^. 

The Luttinger liquid parameters, K p / a , and 
charge/spin velocities u p / a are found by adding 
contributions from harmonic Hamiltonians (|29p and 
£Z£ ter and H?£® from ((311), with the result 

2^(0)-[/(2fc F )cos 2 [ 7F ] , 
up w u F (l H — ), 



u F (l- 



27TU F 

[/(2/c F )cos 2 [ 7F 

27TU F 



K„ 



1 



27TU F 

t/(2fc F )cos[2 7F ] 
27rw F 



(41) 



These expressions are perturbative in small parameters 
U(0)/vf and [7(2fc F )/u F . Note also that physically rea- 
sonable interactions are characterized by U {2k p) < U{0), 
where the equality sign is obtained in the limit of fully 
screened, delta-function like contact interaction between 
electrons. 

Noting that 7F varies from to ir/2 as the ratio 
E S . J A 2 varies from to oo, 



7f = arctan ■ 



2a R k 



rkf 



for E s _ < A 2 
tt/2 for E S . D > A 2 



(42) 



we observe that spin stiffness K a in (|41"j) varies from 
its standard value slightly above 1, K a {j F — + 0) = 
1 + U{2kp)/ {2ttvf), to the value below 1, K a {^F — * 
7r/2) = 1 — C/(2/c F )/(27r-y F ). This unusual behavior, con- 
sequences of which are discussed below, is rooted in the 
spin-orbit-broken spin-rotational invariance of the prob- 
lem, as discussed in the Introduction. 



B. Renormalization Group Analysis 

The fate of the three non-linear terms, Cooper Q37[) . 
backscattering (|38|) and asymmetric (|39[) , are determined 
by renormalization group (RG) analysis. The analysis is 
significantly simplified by expressing the Hamiltonian in 
terms of current operators. To this end, we write the 
Hamiltonian in terms of the right and left spin {Jr, Jl) 
and charge {J R , J£) currents, which obey Kac-Moody 
algebra^. The uniform part of the currents are expressed 
in terms of the chiral right and left moving fermions: the 
charge currents are 

j c R =j2 R t R x, j c l = J2 l » l »> ( 43 ) 

and the spin-currents are 



Jr = 



E Rt-^-R^ 



i Jl — 



E 



L], — - — L h 



(44) 



As an example, the z— component of the right moving 
spin current is defined as J R = {R_R- — rX i?+)/2. Note 
that in the asymptotic limit of a R kp / Az — > the (— , +) 
bands correspond to (t, I) spin bands and we recover the 
canonical definition for the spin current. 

Since the charge part of the Hamiltonian is quadratic, 
(|35[) . and is decoupled from the spin part, it suffices to 
consider the RG flow of the spin part only, H„ = H® + 
+ + H^. In terms of current operators it reads: 



Hi 


= 2%u a J 


H* 


= 


H* 


= ™ a y® 


H C a 


= ™ a ye 



dx\ 



-i28k,FX j— t+ 



j R Jl + h.c.; 



aVa f dx{J R J L +h.C.). 



(45) 



The initial values {£ 
are 



0) of the interaction parameters 



y CT (0) = 2{K a - 1) = U{2k F ) cos[2 7F ]/ 
l£(0) = U{2k F )sm 2 { lF )/TTu <T , 
y*(0)=-*7(2fe F )cos 2 ( 7F )/7ra CT1 
^(0) = C/(2fc F )sin(2 7F )/7ru CT . 



TTVp, 



(46) 
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Note that to first order in U(2k F ) there is no difference 
between v F and u a in denominators of the above expres- 
sions. 

It is convenient to start with formal but useful limit of 
Skp = 0, where H a can be compactly written in terms of 
spin currents 

H a {5k F = 0) = 2ixu a J dx^jyft + J[J[) + 



E 



VaJ a R J a R 



vaWI-J x r jd\. 



Here y x = y& + y% , y y = y^ 
yA = y„- The RG equations for the dimensionless cou- 
plings y a =x, y .z,A are easy to derive with the help of OPE 
(operator product expansion) technique, 



y*, 



(47) 



and 



dy x 
dl 

dy y 

dl 
dy z 
dl 
dyA 
dl 



V-yVz, 

yzyx + y 2 A , 
VvVx, 

y v yA- 



(48) 



Despite complicated appearance, the solution of this sys- 
tem of equations is easy. One finds that y x {€) = y z (l), 
y A {l) = -sin[2-> 'f]v v (£), and y x (£) = cos[2j F ]y y (l). As 
a result, the system is reduced to a single equation 
dy y /d£ = yy, solution of which is standard: y y (£) — 
I/„(0)/(l-y w (0)*)->-l/*for*->oo. 

Thus, in the absence of momentum mismatch Skp be- 
tween the two subbands, all perturbations in (|47[) are 
marginally irrelevant and logarithmically decay to zero. 
This simply reflects rotational SU(2) symmetry of the 
problem in the absence of spin-orbit and Zeeman fields. 



standard system of two KT equations (see, for example, 

Ref EMail). 



dya 
dl 



dye 
dl 



y a yc- 



(49) 



Note that initial values of these couplings are given by 
corresponding solutions of (|48[) evaluated at I = £ z . 

Solution of this system is determined by the integral of 
motion /i 2 = y c — y 2 . and the ratio of the initial couplings 
Vcr(0)/yc(0) = — cos [a]. In terms of these parameters it 
reads^i 

y a (£) = —fi cot [fjd + a] , yc{£) — A*/ sin[/z£ + a]. (50) 

There are three different regimes, illustrated in Figure [3] 

I, strong coupling: 1/V3 > sin[7_p] > 0. Here fi = 
—im, with to > 0, and a = n + i/3, with (3 > 0. In this 
regime both coupling flow to strong coupling, reaching 
pole singularity at Iq = (3/m. 

II, cross-over regime: 1 > sin[7^] > 1/V3- Here \i > 
and < a < n. The flow is still to strong coupling, 
but via an intermediate (cross-over) region (for < a < 
7r/2) where yc(£) initially decreases. Eventually both 
y a _c reach strong coupling, at Iq = (ir — a)//z. 

III, weak coupling: This obtains when fj, — im, m > 0, 
and a = i(3, with (3 > 0. In this situation y a — > —m as 
I — > oo, while yc — > 0. This is critical (Luttinger liquid) 
phase of the spin sector. It is, however, not realized in 
our problem as the requirement — y CT (0) > yc(0) > is 
equivalent to sin^i?] > 1 which is clearly not possible. 

The conclusion is then that Cooper phase is realized 
for arbitrary value of n/2 > > 0, i.e. for arbitrary 
ratio of SO to Zeeman energies, tan[7i?] = 2oLRkpj A z . 
This finding of the Cooper phase, which has the mean- 
ing of the spin-orbit stabilized spin- density- wave (SDW^) 
phase (see next Section), constitutes the main result of 
our work. We have previously discussed the limit of small 
7f, which is physically most transparent, in Ref . l25l. 



yc = -y<r 


/ 


k 2/C VC = Va 






II 










V<r 


FIG. 3: 


RG flow of (l49i 



The full problem, with Skp ^ 0, is solved by neglect- 
ing all momentum non-conserving terms in H a : being 
marginal in the 5kp = limit, such terms become in- 
finitely irrelevant for Skp ^ 0. More carefully, we can 
follow full RG gH) until l z = ln[l/{a Sk F )} = \n[k F /5k F ] 
is reached - beyond this scale oscillating terms average 
to zero. The end result is that we are allowed to dis- 
regard H^' B terms in (|45| . RG equations for remain- 
ing couplings can be obtained from (|48p by nullifying all 
but two, y a and = yc, couplings. This leads to the 



C. The nature of the Cooper ordering 

We now consider the physical meaning of the Cooper 
instability. For simplicity we focus on the regime I of the 
previous Section. Being relevant, the Cooper term (pT7|) 
grows in magnitude and reaches strong coupling limit 
when yci^c) ~ 1 while K a — ► 2 30 . A positive value of gc 
results in 6 a field being pinn ed to one of the semi-classical 
minima 9^ = (m+ \)^-k/2 (to G Z). The energy cost of 
(massive) fluctuations 59^ near these minima represents 
spin gap which can be estimated as 



v F f lF U(2k F ) 
— = v F [ 



irv F 



K a /2(K a -l)) 



(51) 



Here £ = aoe £c is the correlation length, £ ~ 
[Trvp/(U{2k F )j 2 F )] K -/^ K -'- 1 ^. Physical meaning of 
these minima follows from the analysis of spin correla- 
tions. 
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We start with spin density S a = ^\a^ s ,^ s ' /2, which is 
defined with respect to the standard spin basis, s =T, 1- 
We focus on "2fc F " -components of spin, where quotation 
marks are used to remind that large-momentum com- 
ponents of spin density include contributions from both 
k + + k- = 2kF and 2k± processes, see Fig. Q] We find 
(using the gauge n+n- = i) 

S x 

sy 

qz 

15 ' 2k f 



cos[v27riy9p + 2kpx] 
ira 



- sin[V27r6' (T ] 

\ I cos[7_f] cosfv 7 ^^] + sui[7f] cosfV^tyScr + Skpx] 



cos[v27r< ( 5p + 2kpx] 




(52) 



The last line of the above equation is somewhat symbolic, 
with zeros representing exponentially decaying correla- 
tions of the corresponding spin components, S v ' z . Here 
z-component is disordered by strong quantum fluctua- 
tions of dual (fu field, as dictated by [<p, 9] commutation 
relation, see The y-component does not order be- 

cause cosf-s/^^ 1 ] = 0. Thus Cooper order found here in 
fact represents spin-density- wave (SDW X ) order at mo- 
mentum 2k F of the ^-component of spin density, as dis- 
cussed previously in Ref. [25[ Observe that S x ordering 
is of quasi-LRO type as it involves free charge boson, ip p . 
As a result, spin correlations do decay with time and dis- 
tance, but very slowly (S x (x)S x (0)) ~ cos[2kpx] x~ K <' . 

The result (|52| also hints a possibility of truly long- 
range-ordered spin correlations in the insulating state 
of the wire - Heisenberg spin chain. There the charge 
field ip p is pinned by the relevant two-particle Umklapp 
scattering™ which can be mimicked by setting K p — > 
in the spin correlation function above. This is the essence 
of the result to be discussed in Section [VTl below. 

Observe another interesting feature of Eq.152"! S% kp has 
the appearance of rotated by angle component of the 
vector, whereas S% k and S 2k remain unchanged. The 
question that arises is what does S% k rotate into? The 
full answer is provided by considering 2fci?-component of 
the generalized helicity operators 

K kF {x) = -{i/2k F ) Y, ^s<s^A2k F (53) 

= (»/2) Y( R i< S ' L ^- t2kFX - L t< S ' R s' el2kFX ), 

s=T4 

where a = {0, x, y, z}. Note that /iSjfep ( x ) turns into well- 
known staggered dimerization operator e(x) in the 'spin 
chain limit' of the problem, when charge fluctuations dis- 
appear. In terms of the right and left moving fermions 
the staggered (2k f) dimerization is given by 



e=h° 2kF = (i/2) J2( R t L * 
«=T4 



-i2kpx 



L\R f 



).(54) 



Its bosonized form, in terms of <f> a and 9„ fields (|33 
cos[y/2n tfp + 2k px 



n 2k 



7ra 

■ sin[7^] cos[V2T:9 a 



^ cos[7f] cosf-v/^^o- + Skpx] 
(55) 



matches "rotated" S% k in (f52"| exactly. 

Although similar looking, this operator is different 
from ll 2kp v component of the density, described below 
in ([60]) . That one has charge boson cp p appearing un- 
der sine, see ([6"Tj) and (j6"2")l , while both e and S fields are 
proportional to the cosine of it (see also Ref. 41). The 
difference is important. The y— component of spin and 
h operators in the original up and down spin basis can 
be written in a rather compact form, 



ay 

D 2k f 
h° 2kF 



where 



-- cos[7 F ]^ fcF + sin[ 7F ]ft^ F 



\ £ (Rlal^e-^- 
+Ltal v ,R v ,e- i( - k - +k ^ x ) 



(56) 



and 



-Y 

2 ^ 



(RlL v e 



-i2k u x 



LtRvt 



i2k u x- 



(57) 



,(58) 



are, respectively, the spin and h operators in the =F basis. 
Thus, in the limit of 7_f — ► 7r/2, the 2kp component of 
spin along the y— direction in one basis appears as the h° 
operator in the second basis and vice versa. 

Of the remaining staggered operators, h a (a = x, y, z) 
only h y is affected by rotation. The y-component par- 
tially "rotates" into the 2fcp-part of the density opera- 



tor, p 2kF = J2u= T ( R t L " e ~ 
following relation 



i2k v x 



LlR u t 



,i2k v x\ 



h tk F = cos llF}hl kF - sin[ 7F ]p 2fcF /2. 



via the 



(59) 



On the other hand the 2kp component of den- 
sity operator in the original spin basis, P2k F = 
Y< s =-\,i( R l L s e ~ l2kFX + L\R s e akFX ), rotates into the y- 
component of the /i-operator 



P2fc F /2 = cos[7 F ]p 2 fc i! ,/2 + sn\Yf F ]h y 2k 



(60) 



As before, tilde's are used to denote operators in the =F 
basis. The bosonized forms for h\, kp and p 2 k F are as 
follows: 

1 



hn k = sin(v27T0 p + 2kpx) cos(\2~K9 a ) (61) 

and 



P2k f 



■kqq 



sin(v27r0p + 2k F x) cos(V27r0 CT + Sk F x\62) 
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Note that the charge content of h y and density operator 
are the same but the spin parts are different. 

The following relation may be helpful in revealing the 
origins of h a and e fields. In the case of Heisenberg chain 
the staggered dimerization has meaning of the staggered 
energy density, e(x) = e l2kFX S{x) ■ S(x + a ), where 
is the lattice spacing. In the low-energy limit this ex- 
pression turns into e(x) oc (J R (x) + Jl{x)) ■ S2k F {x'), 
where the limit x' — » x must be taken. Short cal- 
culation shows that this leads to Eq. ([54]) above. We 
now observe that, quite similarly to e = h , the helic- 
ity operator h a (with vector index a — x,y, z) may be 
understood as arising from the fusing of the spin cur- 
rent (|44f with the 2/ci?-component of the density field: 
h 2k F ( x ) = ( j r( x ) + Jl( x ))p^ f (x'). Here again x' -> x 
limit is understood. One can check that all relations in- 
volving h a derived above follow from this observation. In 
particular, we note that 



irao 



sin[v27riy9p + 2kpx] sin[v 27r# (T ] 



(63) 



implying that correlation function of this field decays 
with the same exponent (K p ) as that of S% kF discussed 
in the beginning of this Section. It is useful to note that 
helicity disappears in the spin-chain limit of the problem, 
together with the low-energy density fluctuations. 



IV. PERTURB ATIVE APPROACH 

The aim of this section is to show the limit "fp — > can 
be obtained in a straightforward perturbation expansion 
in or. While results of this section parallel conclusions of 
the previous two-subband consideration in Sections IIII1 
MI Al the technical steps involved are somewhat involved 
and are, in our opinion, of interest in its own right. In 
addition, similar perturbative consideration of the impu- 
rity effects later in this work turn out to be very infor- 
mative for understanding the physics. For these reasons 
we choose to present the main steps of the perturbation 
theory in spin-orbit coupling ol r . The calculation starts 
very similar toi^ but concludes with quite different steps. 

The idea is to treat both magnetic field and Rashba 
terms as perturbations to the standard single-channel 
ballistic quantum wire charge and spin sectors of which 
are described by the decoupled Tomonaga-Luttinger 
Hamiltonians (|35[) and (|36|) . The parameters Kp/^jUpu 
of these unperturbed harmonic sectors are given by 
eq. (|4"Tj) but with jp — 0. Spin backscattering term ([55]) 
is in principle present (again with 7^ = 0) but will not 
be required in the subsequent calculation. 

Thus the perturbing terms are, see ([3]), the Zeeman 
term, 



H z = -A z / dx 



(64) 



and the spin-orbit term given by, 

H R = a R J dx *l(x)a* s ,(-i-^)H> s ,(x) (65) 

where \l/ s= |^(a:) right and left movers of unperturbed 
single-channel quantum wire 



* s = R s e lkFX + L s e~ lkFX . 



(66) 



Bosonized expressions for R/L operators parallels that 
in (|22l23l24p where the subband index v = ± should 
be replaced by spin index s =T,|- In terms of charge 
and spin modes introduced in ()33[) , dual pair ip s , 9 S for a 
fcrmion of a given spin projection s is expressed as 



+ sd a )/V2, (67) 



where the following correspondence for the right-hand- 
side of the equations is understood: s =|= +1 and s =1 = 
— 1. It then follows that 



[ Vx {x),B y {x')] = -(1-signOc-O), (68) 



where A = s =f, J, or p, a. 
We then find 



Hz = -A, 
and 



dx (,r R + jd = - 



n/2tF 



dx d x <p a (69) 



H R = 2a R k F J dx (J%(x) - J x L {x)) 
2a R k F n^ l 



na 



J dx cos[y/2~7r(p a } sm[y/2n9 a ], (70) 



where and J£ are the a-th components of chiral (right 
and left) spin-currents (a = x, y, z) defined as 

Jft = R\{x)^fRA x ) , Jl = Ll(x)^fL sl (x). (71) 

Note that H R , being determined by the difference of right 
and left spin currents, is odd under spatial inversion V 
which interchanges right and left movers. 

Finally, we rescale fields 8 a — > 8 a /y/K<j and ip a — > 
V 'K a ip a and account for the Zeeman term ()69|) by a 
position-dependent shift <p a — > ip a + \J K a /(2Tr)A z x/u a 
so that 

r I 2tt 

Hr = g R J dx cos[V 2itK „ ip „ + qox] sin[W — 9 a ], (72) 

2a R k F K a A z 

9R = VlVi, 9o = ■ 

7rao Ua- 

Upto a factor of K a ps 1 which appears here due to the 
rescaling of the bosonic fields above, qo matches with 6k p 
in © for or = 0. 

These standard transformations leave ([72"]) as the only 
perturbation, and correspond to A z ^> a R kp limit of the 
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theory. Observe that charge sector of the theory, (|35|) 
and (|41[) with jp = 0, decouples from the spin sector 
and does not generate any new term. 

The calculation proceeds by expanding partition func- 
tion in powers of gn 

Z = J e~ So {\ - J drH R + 1 J drdr'H 2 + ...), (73) 

where the unperturbed action Sq is particularly simple 

So = Jdxdr ^[{d xVa f + (de a ) 2 } - id rVa d x e a . (74) 



The back-scattering term cos(v87rK^(/3 . + 2Skpx), see 
(I38p . is ignored as it generates terms at higher order in 
perturbation theory (by coupling with H R ). 

Further details of the calculation are summarized in 
Appendix [B The result is (|B12|) fr om where we can 
read off the leading correction to the spin Hamiltonian 



/ a R k F s 


,2 u{2k F ) f 




1 (na ) 2 K a J 



- f dxcosJ^e*}. (75) 

cr J V 



This compares well with our previous result (|37p (remem- 
ber that we rescaled 8 as 8 — ► 8/^K above) in the limit 
7f <C 1, where the described calculation is applicable. 

It is worth pointing out that perturbatively gener- 
ated Cooper term, instead of being proportional to 
(aftkp/vp) 2 as one would naively expect from a straight 
forward perturbation expansion of (|72|) , acquires a non- 
trivial dependence on both the back-scattering ampli- 
tude, U(2kp)/ (2itvf), (via K — l/K combination in 
(jBlip ) and the Zeeman term A z (via qo dependence of 
(jBlip ) and is proportional to (a^kp / A z ) 2 U(2kp)/vp. A 
second thing to notice is the crucial role of Klein-terms 
(see (1B12[) ) in generating the correct (positive) sign in 
ea.([75|). 



A. Arbitrary angle between SO and magnetic field 

Perturbative approach is convenient for analyzing an- 
gular stability of the Cooper phase. Let us consider the 
situation when the angle between the magnetic field and 
spin-orbit directions is ir/2 — (3, so that (3 = corresponds 
to the orthogonal orientation studied so far in this paper. 
It is convenient to keep z as the direction of magnetic 
field, in which case the spin-orbital term changes from 
ajipa x to ajip(a x cos[/3] +<r z sin[/3]). Following steps that 

led to (|69I70|) and ([72]) we obtain that Hr,p 



f I 27T 

H R = 9r J dxcos[y/2irK rT ip cr + q x] sin[W j^-^a], 
Hp ] = -<J^a R k F sm[(3] J dx d x 8 a , (76) 

where first (second) equation represents contribution 
from a x (cr 2 ) matrix. The coupling constant of cosine 



term is only slightly modified in comparison with (|72p . 

Sr^ = 9R cos2 [P}- Hr^ represents the main new feature 
of the non-orthogonal situation. It is naturally absorbed 

(x) 

into H R by a simple position-dependent shift of field 8 a 
which results in 

Hr,{3 = J dx cos{^2iTK (J tp a + q x] 

2aRkp sin[/3] 



x sin[W — d a + Sx] ; S = 
K, 



llT^vp 



(77) 



Note that both ip and 8 fields acquire oscillations now. 
The consequence of this is that the generated Cooper 
term (compare with ([75]) ) 




(78) 



does not conserve momentum for (3^0. This important 
result is evident in the single particle spectrum, which 
now reads (/3 ^ 0) 



2m 



± V(<*Rk) 2 + (A./2) 2 - akA z sin/3, (79) 



and is to be compared with ([4]), which describes /3 = 
situation. The spectrum (|79p acquires an asymmetry 
about the energy axis as seen in Fig. [U The top (+) 
subband shifts towards the right and the bottom (— ) one 
towards the left. As a result, the pair-tunneling docs not 
conserve momentum for (3 ^ 0, as is seen in ([78]) above. 




FIG. 4: Occupied subbands e± for the case of non-orthogonal 
spin-orbital and magnetic field axes. Arrows illustrate spin 
polarization, as in Fig. [T] Filled dots indicate location of 
the center-of-mass for (+) (red) and (— ) (blue) subbands. 
Dashed lines show (±) subbands of Fig. Q] corresponding to 
the orthogonal orientation, f3 = 0, for comparison. 

The momentum-mismatch S de-stabilizes the Cooper 
order and eventually, for some critical f3 c , destroys it 
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completely. The critical angle can be easily estimated by 
comparing two spatial scales: £c> which describes per- 
fect SDW order, and 1/5, which represents the scale on 
which momentum-nonconservation becomes pronounced. 
Estimating £c ~ eia , where gn(i c ) ~ 1, we obtain that 
the Cooper order is destroyed when 

U(2k F )(a R k F cos[/3]/A z ) 2 » a R k F sin[/3]. (80) 

Observe that this happens already for small angles, we 
can simplify the expression for the critical angle 

C oc U ff a £* « 1. (81) 

This angular sensitivity of the found Cooper order to 
mutual orientation of the magnetic and spin-orbital di- 
rections can be used as an experimental probe to dif- 
ferentiate between it and other, spin-orbit independent, 
many-body instabilities of interacting quantum wire. 

V. ISOLATED IMPURITY 

SDWa; state manifests itself not only in spin correla- 
tions (|52p . It turns out that its response to a weak poten- 
tial scattering (impurity) is rather non-trivial. We con- 
sider here a weak delta- function impurity V(x) = Vq6(x), 
with strength Vq, located at the origin. The condition 
Vq <C A c means that impurity can be considered as a 
weak perturbation to the established SDW phase. Since 
the latter is most robust in the limit of -f F <C 1, this 
is the limit we consider in this Section. (The limit of 
strong impurity, Vq 3> A c , is rather standard: impurity 
destroys the SDW and the wire flows into an insulator at 
low energies 3 -!.) 

The interaction of electrons with an impurity potential 
V(x) is given by 

V = f dx V{x) *l(x)* a (x). (82) 
■* s=T4 

As usual, it is the backscattering {p2k F ) part of the den- 
sity that has to be considered, p2k F = T1 S {R\L S + hue). 
Working in the two-subband basis of eigenstates v = ± of 
the non-interacting Hamiltonian which includes Rashba 
spin-orbit term and the Zeeman term, see (|14p. we derive 
for the 2fc^-component of the density p(x) at the origin 

P2U F =P\T + PT;, (83) 
where the inira-subband part is 

P2kT = cos[ 7F ](LlR + + LlR- + h.c), (84) 

while the novel inter-subband contribution is present due 
to the non-orhogonality of spin states in (±) subbands 
discussed in Section [H] 

PiV = sHlF]{LlR + - L\R_ + h.c). (85) 



Alternatively, one can think of these two contributions 
as representing p and h v contributions in (|60[) . In terms 
of bosonic fields the density reads 

— 2 sin[v / 27n ) s J / i — 
P2k F (x = 0) = ^ cos 7f cos V27r¥V 

7TO0 V 

- sin[7 F ] cos[V27^ CT ]) • (86) 

We now observe that setting a — > 9f = (m ^-^)-JtF/2, 
as appropriate for the Cooper phase (Section [ill C[) , nul- 
lifies the backscattering component of the density (|86|) . 
The first (intra-subband) term gets killed by diverging 
fluctuations of tp a , dual to 6 a . Intriguingly, the sec- 
ond (inter-subband) contribution is also zero because 
cos[V2^9f] = ±cos[tt/2] = 0. 

This argument can be made more precise by fol- 
lowing calculations described in Refs. SHHi. We 
parametrize a = Of + SO and expand the relevant co- 
sine term in Hamiltonian ([57)1 to second order in fluc- 
tuations 69. One obtains a massive term oc A c (69) 2 
in the Hamiltonian, which causes exponential decay in 
correlation functions of the dual ip„ field. In par- 
ticular (cos v / 27t</' ( t(0, t) cos y/2nip a (Q, t')) will decay as 
exp[— A c |t — t'\]. Thus for an incoming particle with 
energy u> <C A c , the inter-subband scattering channel is 
absent. 

Substituting 9„ = Of + 69 in the second term in ([5S|) 
converts it into ^ F sm[\Z2n(p p } sin[\/2n 60]. Correlations 
of sin[v / 27r<5^] are also short-ranged 

(sm[V2n56(T)] sin[V2^6>(-r')]) « x 

v F k F 

sinh[X (A c |r-r'|)], (87) 

where Kq(x) ~ e~ x / ' \Jx (for x 3> 1) is the modi- 
fied Bessel function (see Refl29[). The two exponen- 
tially decaying contibutions add up (in second order 
perturbation theory in impurity strength Vo) to pro- 
duce an effective two-particle backscattering potential oc 
(Iq 2 /A c ) cos[v / 87T( ( 5p(0, t)]. This generated two-particle 
impurity backscattering term, however, is relevant only 
for strongly repulsive interactions, K p < 1/2. We are 
thus left with irrelevant impurity potential for as long as 
Vq -C A c is justified and for not too strong repulsion, 
1/2 < K p < 1: SDWa, state is not sensitive to weak 
disorder! 

The situation is similar to that in recently proposed 
edge states in quantum spin Hall system 3 -^. There, gap- 
less spin-up and spin-down excitations propagate in op- 
posite directions along the edge, which forbids single- 
particle backscattering. Interacting electrons, however, 
can backscatter off the impurity in pair o 38 ' 39 . 

Our conclusion p2k F ~ * 0j ea. (|86[) . rests on somewhat 
technical condition cos[V2ttO^] — and deserves a bet- 
ter understanding. This is provided by the perturbative 
calculation in Appendix [DJ The idea of the calculation 
is similar to that in Section [TVl treat both impurity and 
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spin-orbit terms as perturbation and generate spin-orbit- 
related corrections to backscattering potential (propor- 
tional to 7f in (|86|) perturbatively. In this way we can 
be certain that all symmetry-allowed contributions are 
accounted for. This argument also makes it clear that 
the generated terms, being produced by the SOI, have 
to be odd under spatial inversion V . It is this symmetry 
that guaranties that sin[v / 2~7r6' -] can not be generated in 
the process. This instructive calculation, carried in Ap- 
pendix 02 indeed supports the conclusions of the more 
formal two-subband approach described in this Section. 

We conclude this Section by pointing out experimen- 
tal consequences of our findings. Suppression of single- 
particle backscattering off weak impurity under the out- 
lined conditions implies an unusual negative magnetore- 
sistance in one-dimensional wire. Indeed, the conduc- 
tance of the wire with such an impurity should remain 
at the perfect G w i ro = 2e 2 /h value for as long as ap- 
plied magnetic field is directed perpendicular to the spin- 
orbital (a x here) axis. By either turning the magnetic 
field off, or simply rotating the sample (so that mo- 
mentum mismatch between the two subbands destabi- 
lizes SDW X order) , one should observe that conductance 
plateau deteriorates. It will be destroyed completely in 
the zero-temperature limit. 



VI. HEISENBERG ANTIFERROMAGNET 
WITH DM INTERACTION AND MAGNETIC 
FIELD 

In this section we will consider a closely related prob- 
lem: the effect of uniform and staggered Dzyaloshisnkii- 
Moriya (DM) interaction on a spin S = 1/2 Heisen- 
berg antiferromagnet (HAFM) in the presence of a mag- 
netic field perpendicular to the DM vector. We start by 
demonstrating the equivalence between the quantum wire 
problem, analyzed in the previous Sections, with that of 
Heisenberg spin chain subject to an additional uniform 
DM interaction. We then present a novel chiral rotation 
argument to show that results of the two-subband ap- 
proach in Section IIIII can be obtained from a straightfor- 
ward combination of two independent rotations for right- 
and left-movers. We also discuss the difference between 
the uniform and the staggered DM interaction, analyzed 
previously in Ref. [27]. 



A. Weak uniform DM interaction and strong 
Magnetic Field 



Zumino-Novikov-Witten model the basic ingredients of 
which are the uniform SU(2) left Jl(x) and right Jr(x) 
spin currents, defined via fermions in (|7ip . and the stag- 
gered magnetization N(x)^. Another important compo- 
nent is provided by staggered energy density (dimeriza- 
tion) , see eq. flMJ) , but this does not develop any order in 
the presence of external magnetic field. These operators 
will be used to represent the continuum limit of the spin, 



a (j L {x) + J R {x) + (-l) x / a °N(xj), (89) 



where oq is the lattice spacing and continuous space vari 
able is introduced via x = ja$. Continuum limit of 



H 



2nv 



dx Jr ■ Jr + Jl ■ Jl 



- I dx (d x <p) 2 + (d x ^ 2 



(90) 
(91) 



The first line constitutes nonabelian spin current formu- 
lation of the problem, which is convenient for analyzing 
spin-rotation invariant (SU(2)) problems, while the sec- 
ond makes connection with familiar abelian bosonization 
result (|36| : identifications {ip a ,O a ) — > (<p,0), u„ — > v = 
irJao/2 and hp — » 7r/(2a ) finalize the connection. The 
marginal backscattering term accounts for residual inter- 
action between spin excitations, 



-ffbs = 



dxJR ■ Jl- 



(92) 



Its coupling constant c/bs > is known from extensive 
numerical studies of the Heisenberg chain^ : <7bs = 0.23 x 
2irv. This term is responsible for the fact that initial 
value of the Luttingcr parameter K a = 1 + gb s /2Trv is 
greater than 1. 

There are no low energy charge fluctuations as the 
charge is "locked" to the lattice by relevant Umklapp pro- 
cesses. This allows one to replace cosine of charged boson 
ip p in (|52"|) by its expectation value A: A = (cos[\/2iripp]). 
The first line of (|52[) then establishes abelian representa- 
tion of the staggered magnetization 




A(-l) 3 
na 



- sin[V27r6>] 
cos[V2tt0] 

— sin[\/27riy9] 



Magnetic field is introduced via Zeeman term 



H 



J2 WeBS* - —A z J dx(J R + Jl), 



(93) 



(94) 



The Hamiltonian of an isotropic HAFM spin chain is 
given by 

-ffhcis — J ^ ^ [Sj Sj+i + SjSj^-^ + SjSj^i]. (88) 



In the continuum limit, the above Hamiltonian is most 
conveniently described in terms of the SU(2)i Wess- 



which is identical to (|69|) . The small upper index {z 
here) indicates the axis in spin space. Finally, the 
Dzyaloshinskii-Moriya Hamiltonian describes asymmet- 
ric and odd under spatial inversion interaction between 
spins 



^DM — D ■ S-j X S 



(95) 
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Vector D = Dx fixes direction of spin anisotropy, which 
we choose to be along spin-i direction, similar to the 
spin-orbit direction choice in Q. 

Using (|89p we see that continuum limit of the uniform 
DM Hamiltonian requires the knowledge of the following 
operators 

Oj(x,x') = J y (x)J z (x') - J z {x)J y (x r ), 
N (x,x') = N y {x)N z {x') - N z (x)N v (x'), (96) 

where x' = x+ao- The first of these follows from the well- 
known OPE for spin currents, see for example eq. (25,26) 
of Ref. Hi] and set z — z' = ia : 

0j (x,x>) = —(r R (x)-r L (x)). (97) 

ira 

The second requires more work and is calculated using 
the definition (|53J) and bosonic OPEs [[EI]) 

N (x,x') = — (J%(x)- Jl{x)). (98) 
7ra 

As a result 

#dm =Djdx (J% - Jf) , D = ^(1 + 2A 2 ). (99) 

Comparing Q with (JTOJl and ([M]) with |BS]) we observe 
that the problem at hand is identical to that analyzed in 
Section ITVl We immediately conclude that the 2kp = n 
(staggered) components of spins acquire long range or- 
der along the direction of D as soon as magnetic field is 
turned on along the orthogonal direction. 

(N x ) = const (-If, (N y ' z ) = 0. (100) 

It is useful to compare this equation with ((52)) . Note that 
the effect is controlled by the backscattering amplitude 
g hs = 2{K a - 1). 

Except for the brief remark in section VIII of Refill, 
the possibility of the long-range order in the magnetized 
spin chain with asymmetric uniform DM interaction has 
not been discussed in the literature, to the best of our 
knowledge. We also would like to note technical similari- 
ties of our problem with that of an anisotropic Heisenberg 
chain in a transverse field, considered in Refl43l 

Experimentally the field-induced long-ranged mag- 
netic order would be probably easiest to observe via ther- 
modynamic measurements. Similar to the case of copper 
benzoate^, the order should show up via an exponential 
suppression of the specific heat. The magnitude of this 
suppression, which is controlled by the energy gap A c 
(f5Tj) , should be a sensitive function of the magnetic field 
orientation, reaching maximum when the field is orthog- 
onal to the DM axis as discussed above. 



B. Chiral rotation and tilted magnetic field 

We now present an elegant argument, borrowed from 
the technically closely related study of quantum kagome 



antiferromagnet , exposing the nature of the Cooper or- 
der to the fullest. Consider the situation, discussed in 
Section IIV Al of the tilted magnetic field making an an- 
gle 7r/2 — (3 with the DM axis direction (which we keep 
at x). The Hamiltonian describing this arrangement is 
given by H = H + f/bs + V, where V includes Zeeman 
and DM fields 

V = J dx(d R J x R - d L J'i - h x {J R + J£)). (101) 
The parameters are 

d R — D — h%, dp = D + /i2 ; hi = hcos/3, hi = /isin/3. 

(102) 

We now rotate right (left) spin currents Jr (Jl) about 
y axis by angle 9 R (0l) such that after the rotation the 
V-teim becomes 




The relation between old (J) and new (M) currents is 
simple 

Jr = K(9 R )M R , J l = K(e L )M L , (104) 

here 1Z is the rotation matrix 

/ cos[0] -sin [9] \ 
Tl{9) =01 . (105) 
\ sin[0] cos [61] / 

The rotation angles are given by tan[0^] = d R /h\ and 
tan[0£,] = —dh/hi. The key reason behind these rota- 
tions is the observation that unperturbed Hamiltonian 
Hq, being the sum of commuting right and left terms, 
is invariant under the rotations. But the backscattering 
-f/bs is not, and transforms into 

H hs = -g hs J di(cos[ 7 ](M£M£ + M R M[) 

+ sin[ 7 ](M£M£ - M%Mt) + M y R Mf) (106) 

where 7 = 9 R — 9l- The reason for this transformation 
is of course that right and left currents are rotated in 
opposite directions (and by different amounts). Observe 
that (|106[) matches the interaction part of (|4T[) . 

The nice thing about the rotation is that now both 
SO and Zeeman fields can be taken into account by sim- 
ple linear transformations of spin bosons. Indeed, using 
abelian bosonization we observe that 

V = - [ dx(^Ld xVa + ^Ld x 9„), (107) 

where 

t v - (y/cP L + hf + yl# R + h\)/2v, 

te = (y/cP L + hl - ^Jd 2 R + hl)/2v. (108) 
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These linear terms are removed by shifts 

<P<7 -> fa 



/2tt 
tex 

V2n' 



(109) 



The price is that transverse components M R J L acquire 
oscillating position-dependent factors 



M+e^ +tt) ". (110) 



The major consequence of this is that every term in -ffbs, 
with a single exception of MjjM£ one, picks up oscillating 
factor 



Hbs — —gbs 
cos[7] — 1 



4 

cos [7] 



dx( cos [7] M R M Z L + 



(M+M+ e* 2 *^ + h.c.) 



1 



, (M /r V/ ; ;e ' 



h.c.) 



(Ill) 



+ ^n^ M z M + e -i(t v -t e )x _ M z R M+e i( -^ +tl > )x + h.c.)) . 

This equation contains all arrangements that we have 
discussed in this paper. Setting tg = corresponds to 
orthogonal ((3 = 0) orientation of spin-orbit (DM) and 
magnetic field axis. In that limit we recover equations 
(|45|) (J there corresponds to M here). Allowing for j3 > 
leads us to Section IIV A[ results of which represent 
abelian version of the discussion here. As shown there, 
the SDW X phase is stable in a finite angular interval near 
,3 = 0. 

It is worth noting that (|103[) and (|107[) imply finite 
expectation values of M R ^ L currents: 

(M Z R/L ) = y/<P R/L + hl/(4nv). (112) 

By virtue of (|105[) this implies finite expectation values 
of the uniform magnetization along z and x axis, J R 



/x 



, and uniform spin current, J R X 
two axes. That is, 

{Jr-JD 

(Jr + Jl 



J? , along these 



Z7T 



cos[ 7 ]^, 

trh 



^[7]^,(Jfl + J£)=cos[ 7 ]^.(113) 



Equilibrium values of spin current and magnetization 
along the y— direction vanish. These relations comple- 
ment discussion of the relations between staggered (2k p) 
components of various fields in Section IIII CI 



C. Strong uniform DM interaction and weak 
Magnetic Field 

It is instructive to consider another 'tricky' limit of 
the problem: strong uniform DM interaction and a weak 



magnetic field, D ^> A z . The idea here is to account 
for the DM term exactly on the lattice level, see (II 14[) . 
and treat the Zeeman term as a small perturbation to 
the obtained Hamiltonian Q116p . 

The two axis, DM and Zeeman, are still orthogonal but 
it is more convenient now to let D point along the z-axis 
in spin space, while magnetic field is pointing along x. 
Following Ref. [42] we perform unitary transformation of 
the lattice Hamiltonian to absorb the uniform DM term 
exactly, 



st ste iaj , s: 



, ^ and 5| - 5;, (114) 



where a = arctan(D/J) ps D/J for D <C J. The trans- 
formed Hamiltonian is that of the XXZ spin chain with 
weak anisotropy 

h = + s]S] +1 + is*s* +1 ], (115) 
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where J = \/ J 2 + D 2 . Magnetic field term, however, 
acquires position dependence under this transformation 



H 



(x) 



■ 5^(5/ e iai + Sje- ia] ). (116) 



Bosonizing it according to the rules described in Sec. IIII Al 
and IIVI we obtain 



H 



(x) 



7TO 



dx sin[v27rif(^] cos 




a:Ep-17) 



where Luttinger parameter of the XXZ chain (|115[) is 
given by R- 1 = l-arccos[J/V^ 2 + D 2 ]/n « l-D/(nJ). 
The shift of 9 field by ax can be easily understood from 
equation (|9"9"|) , adapted for the DM axis along i-direction 

in spin space: ifj,^ = D J dx(J R — J£) cx D J dx d x 9 
which indeed can be accounted for by shifting 9. 

Being of nonzero conformal spin, the Zeeman term 
(| 11 7(1 generates under RG two new terms of importance 
to us: H v oc (A 2 D/ J) cos[y/8Tr Kip], of scaling dimension 

2K = 2 + 2D/(irJ) > 2, and H e = G e J dxcos[^9 + 

2ax] with scaling dimension 2/K = 2 — 2D/(xJ) < 2. 
The position-dependent phase 2ax, nonetheless, makes 
it irrelevant. 

Indeed, the coupling constant can be estimated, fol- 
lowing momentum-shell RG of Appendix |C| as Ga oc 
(A z /v) 2 (K - 1/K) 2 » (A Z D/J) 2 , see (JC7| and dCT2l) . 
Neglecting the oscillating phase for the moment, we 
can easily estimate the correlation length correspond- 
ing to the Cooper order as £ c = exp[£ c ], where 
G e (Q = G e (0)exp[(2 - 2/K)Q w 1, so that i c = 
(7rJ/D)bx[J/(A z D)]. 

This length is to be compared with £ a = 1/a w J/D, 
which is the length on which the cosine in Hg changes 
sign. We observe that 



ttJ 
~~D 



In 



J 



A Z D 



(118) 
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which implies that fast oscillations make it impossible 
for the Hg to reach strong coupling. Hence no order can 
come from this term. 

The other term, H v , appears to be irrelevant and it is 
tempting to disregard it altogether. One however must 
be careful and recall discussion in Section IIII Bl Par- 
ticularly, observe that H v may fall in region II (with 
y a {fi) < 0). Whether or not this happens is in princi- 
ple a subject of precise calculation which however can 
be avoided here. Indeed analysis of Section IVIB1 taken 
together with results of Sec. IIII Bl shows that Cooper 
instability develops for arbitrary 7 as long as t$ = 0. 

We thus conclude that H v must belong to region II 
in Figure [3J even though this conclusion is not at all 
obvious from the abelian bosonization analysis sketched 
above. We chose to present this subtle point in order to 
demonstrate the power of non-abelian current formula- 
tion in Section IVlBI 

Note finally that H v still describes Cooper order even 
though it is written in terms of ip field. The reason for this 
illustrates another shortcoming of abelian formulation. 
Unitary transformation (I114j) is the lattice version of ro- 
tation of right (left) currents by 6r — ir/2 {9l = — 7r/2) 
in the previous Section. (In this limit A z = and mag- 
netic field is to be added later as a perturbation.) Such 
a rotation corresponds to 7 = n which replaces standard 
backscattering (M~^M£ + h.c.) with (m£m£ + h.c). In 
terms of abelian bosonization this corresponds to inter- 
change cos \/8ttip <-* cos \ZStt9. Moreover, this same nota- 
tion changes the sign of MSM£ term, which corresponds 
to interchange K <-» 1/K, which completes the duality 
mapping. Observe again the ease with which current for- 
mulation of Section IVI Bl leads to conclusion. 




FIG. 5: Band structure in the absence of magnetic field. Note 
the absence of the gap between the subbands at k — 0. 

In simpler terms, accounting for DM (or, equivalently, 
spin-orbit) term only amounts to working in the basis 
of eigenstates of operator a x . These eigenstates are de- 
scribed simply by two horizontally shifted parabolas in 
Figure [5l Denoting the right (left) subbands by index 

1 (I 1 ) = I ->•)) and 2 (I 2 ) = I «-)) we immediately ob- 
serve that Cooper scattering discussed previously cor- 
responds to the process like R\LiL\r 2 . (So that ini- 
tial pair L1R2 scatters to final R1L2 state. Note that 
the other possibility, L\ — > L2 and R2 — > R\ is for- 
bidden by orthogonality of single-particle basis states, 
(1|2) = (— > I <— ) = 0.) Bosonizing this in a standard 
way we find that (recall R\Li ~ exp[—iV^br(0ii? + </>il)]) 
R[L 1 L\r 2 ~ exp[— i-\Z4jr(<pi - tp 2 )] = exp[—iy/&7rip„]. 



Thus indeed the Cooper process is written in terms of 
ip a field in the new basis. 

Observe that this conclusion is special to the limit of 
no magnetic field, A z = 0. As soon as the field is on, 
A 2 7^ 0, the band structure changes and horizontally 
shifted subbands (1,2) turn into vertically shifted pair 
(— ,+) used in this paper (Fig. [T|): arbitrary small B 
leads to the gap A z at the point k = 0, see ((U). We thus 
conclude once again that A z = is a rather singular 
limit, not much suited for the consideration of general 
situation with both spin-orbit and Zeeman fields present. 



D. HAFM chain with staggered DM interaction 
and magnetic field 

For the sake of completeness and uniformity of the pre- 
sentation, we re-visit here the well-studied case of in- 
terplay between staggered DM interaction and Zeeman 
magnetic field, described in Refs. I27I and[45l 

Staggered DM interaction along ai-axis is described by 

H^m = DY,(^y{S v 3 S z ]+1 - S z S y +1 ). (119) 

i 

Its continuum limit requires the knowledge of 

Oj N (x,x') = N y (x)J z (x') + J z (x)N y (x') 

-J v {x)N z {x') - N z (x)J y (x'),(120) 

where, as before in we set x' = x + o,q. It is 

an easy calculation, using eq.26 of Refl4ll. to find that 
in the absence of magnetic field Ojn = identically 
The situation changes when magnetic field (along z- 
axis) is present. Using bosonization, magnetic field is 
accounted for (see (IT2"j) and line above it) by the shift 
y/2TrK(f — * \I2-kKlp + qox. Moreover, the main effect of 
the field here is contained in qo = A z /w and we can keep 
K = 1 in all calculations below. Observing that a shift in 
ip implies equal shifts in 4>ril fields (so that 9 does not 
change) and working backward through ([52"]) . (fTTj) and 
(122]) we conclude, following RefH3, that in the presence 
of the magnetic field, spin excitations along the field (z- 
axis here) and transverse to it (x, y axis) have minima at 
different momenta. Namely, while J z is still centered at 
5 = 0, transverse components of spin current J x ' y acquire 
minima at q — ±qq — ±27rm, where m is the magneti- 
zation. In addition, N z shifts from n to q — n ± 27rm 
while staggered transverse components N x ' y remain at 
2k p = tt point. As a result, the product N V J Z (first line 
in (|120[) ) retains its zero- field structure and continues to 
remain at zero. At the same time the other combination, 
J y N z (second line in ^20]) ). splits into slow (e^oO*-*')) 
and fast (e l <ia( x + x )) oscillating pieces which do not can- 
cel each other anymore. The remaining calculation is 
most conveniently performed using fermionic representa- 
tions of spin currents (I7ip and longitudinal magnetization 
N z = |[-R}L T - L\R L ]e- iqox + h.c. Fusing right (left) 
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movers of like spin using (|D4[) we obtain, for example, 

(J+(*)+J+(*))JVV) = - Si "L^~n )] (4 L l+ L l^)- 

(121) 

Bosonizing this expression (note that it is not sensitive to 
the sign of the coordinate difference (x — x')) we finally 
obtain 

Oj N (x, x') = -t^- a»[V^6(x)]. (122) 

The same result can be obtained, after somewhat longer 
calculation, using bosonized forms of spin currents and 
magnetization from the very beginning. The continuum 
limit of staggered DM term then follows 

H sDM = / dx cos[V2^(a;)]. (123) 

This is a highly relevant operator (scaling dimension 
1/2), the coupling constant G s -dm of which grows as 
G s -dm{£) — G s _dm(0) exp[3£/2]. This, and the depen- 
dence of its initial value on the combination DA Z (which 
enters via dependence on qo), leads to the energy gap 
in the system that scales as (DA Z ) 2 ^ 3 , exactly as Refl27l 
found originally. Note finally that (|123[) and (|52")l implies 
that the spins order (in a staggered way) along y-axis, 
orthogonal to both DM and magnetic field directions. 

Note that although we have treated staggered DM 
term (|1 19|) as a perturbation to the spin chain subject 
to magnetic field (which comes in via the momentum 
go here), we obtained the same strongly relevant result 
(ETJgJl as the authors of Ref. did. But Ref. arrived 
at it from a different limit: the authors used staggered 
version of (|114p to account for DM piece (|119|) exactly 
and then added magnetic field as a perturbation. The 
outcome (|123p is obtained in both cases simply because 
it is more relevant (dimension 1/2) than the magnetic 
field (Zeeman) term (dimension 1): no matter how small 
Gs-dm initially is, it controls the physics. The 7V y -order 
is present for any ratio of G s _dm/A z ratio. 

VII. IMPLICATIONS FOR ESR EXPERIMENTS 

In recent years, electron spin resonance (ESR) tech- 
nique has become a leading experimental candidate for 
probing anisotropic terms in the spin chain. Recent the- 
oretical impetus to this field has been provided in an 
important work by Oshikawa and Afflec k 46 ' 47 who dis- 
cussed limitations of earlier theoretical wor k 48 i 49 i 50 and 
improved and extended upon them by using powerful 
modern theoretical techniques. The emphasis in^ 7 - has 
been to study the role of anisotropic terms, in particular 
the staggered DM interaction, see Section IVID1 and ex- 
change anisotropy terms, in modifying the resonance po- 
sition and the line width. However, ESR in a spin chain 
with uniform DM interaction was not considered. We will 



follow a closely related work by De Martino et a L 51 ' 52 , 
who considered ESR for carbon nanotubes with spin- 
orbit coupling, to investigate the modifications brought 
upon by the uniform DM term on the ESR spectra. 

We consider Faraday configuration in which the static 
magnetic field and oscillating field are orthogonal to each 
other. The oscillating driving field (of frequency to) is 
in the microwave frequency regime and for all practical 
purposes the spatial modulation of the field can be ig- 
nored. The ESR intensity is given by the transverse (to 
the static magnetic field) spin structure factor at q = 
and frequency u>, 

J(w) = J dtdx e Mt J2 (Jr(x,t)J~,(0,0)). (124) 

r,r'=R/L 

As before, the Hamiltonian, H = Ho+Hb s +V, consists of 
the free part ([50]) . the back scattering part (fi?2")) . and the 
DM and Zeeman terms contained in the V term (1101 1) . 
The angle between DM and magnetic field directions is 
w/2-/3, see Section ECS 

Our aim here is to present the basic picture of ESR 
response for the case of the spin chain with uniform DM 
interaction (equivalently, quantum wire with spin-orbit 
interaction). For this zeroth order description we omit 
the backscattering interaction ([52")) between the spin cur- 
rents altogether. Under this drastic approximation right 
J R and left Jl spin currents are decoupled. This implies 
the intensity (|124|) is the sum of right and left contribu- 
tions, 7(w) = Ir{lo) + 

To account for the simultaneous presence of DM and 
Zeeman terms, we now rotate the right and left currents 
about the y— axis as described in Section IVI B[ see in 
particular (|104|) and (|105j) . As the backscattering p06p 
is neglected, the full Hamiltonian is given by 

H= fdx ^(M« R M a R + MlMl) 

a—x,y,z 

-X R M R - X l M z l , (125) 



where X R/L = Jd 2 R / L + h\ = v(t v T U)- 

We now focus on the contribution of the right spin 
currents, 

Ir(u) = J dtdx e iut {J+(x, t), J R (0,0)), (126) 

which, in terms of the rotated currents, reads 

I R (w) = J dtdx e lu ' t (^cos 2 e R (M^(x,t)Mf l (0)) (127) 
+{M y R (x, t)M y R (0)) + sin 2 R (M R (x, t)Mjj(O))) . 

The cross-terms of the kind {M%{x, t)M R (0)) = due to 
the absence of coupling between x and z components in 
the Hamiltonian (|125|) . Further, switching to M R com- 
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binations of the currents, the intensity becomes 



VIII. CONCLUSIONS 



+ 



+ 



dtdx e ZbJt ( sin 2 6 R {M a (x, t)M z R (Q)) 
(cos 2 R - 1) 



^cos 2 6 R + 1) 



-(Af+(x,i)A/+(0) + h.c.Xl28) 



(M+(x,t)M R (0)+h.c.) 



The second line in this expression contributes zero 
as it requires anomalous averages of the kind 
(* fl , T (a;,i)*i ?iT (0,0)) which are absent in (|125|) . We now 
absorb "right" magnetic field Xr via the shift of right bo- 
son so that -> M^e" tXRx/v (compare with (fTTUj) ). 
It is worth noting that by rotating the spin currents we 
have mapped the problem of the spin chain with uniform 
DM and magnetic fields to that of the chain in field A_r 
(Al) for its right (left) moving components. This allows 
to borrow results of Ref. |4?] and conclude that the first 
line in (| 1 29|) does not contribute to In while the last line 
gives 



(cos 2 0r+1] 



■lu6{u - \Xr\). 



(129) 



Clearly the contribution of the left-moving sector is ob- 
tained by replacing R — > L in the expression above. Thus 
the ESR signal consists of two sharp lines, as previously 
discussed for the case of carbon nanotube in Refs. [5lll52l . 

= (cos 2 6 r + 1)\ r 5(lu- \X r \) 

+ {cos 2 9 l + 1)X l 5{lo~ \X l \). (130) 

The distance between the lines is Al — Xr = 2vtg. The 
relative strength of the two lines is 



Ir(Xr) 
Il(Xl) 



hi 



4 



hf 



(131) 



The ratio is always 1 for (3 = (orthogonal orientation of 
DM and magnetic field axes), see (|102p . Note that this 
is exactly the configuration in which SDW order devel- 
ops, as described in detail in Section UlIBI The ordering 
is driven by the Cooper term, in (|45|) (equivalently, 
the term proportional to (cos [7] — 1) in (|112p ). which 
comes from the backscattering process as is made clear 
by the discussion in Section IVIBI The result (|130[) . ob- 
tained by neglecting backscattering, is clearly not appli- 
cable at low temperature where the SDW (Cooper) in- 
stability develops. Well below the ordering temperature 
the system is described by the sine- Gordon model exci- 
tations of which are massive kinks^. But even above the 
ordering temperature in the disordered phase (which, for 
(3 > extends all way down to zero temperature, see 
Section TlV Ap we expect the backscattering to affect the 
ESR signal. Whether or not it would lead to the finite 
linewidth of the two lines in (| 130[) we do not understand 
yet and leave this interesting question (see Ref. |47] for de- 
tailed discussion of some technical subtleties) for future 
studies. 



Spin-orbital interactions result in reduction of spin- 
rotational symmetry from SU(2) to J/(l) in one- 
dimensional quantum wires and spin chains. This reduc- 
tion, however, is not sufficient to change the critical (Lut- 
tinger liquid) nature of the one-dimensional interacting 
fcrmions. The situation changes dramatically once exter- 
nal magnetic field is applied, as we have shown in this pa- 
per. Most interesting situation occurs when the applied 
field is oriented along the axis orthogonal to the spin- 
orbital (or, Dzyaloshinskii-Moriya, in case of spin chain) 
axis of the wire. The resulting combination of two non- 
commuting perturbations, taken together with electron- 
electron interactions, leads to a novel spin-density- wave 
order in the direction of the spin-orbital axis. 

The physics of this order is elegantly described in terms 
of spin-non-conserving (Cooper) pair tunneling processes 
between Zeeman-split electron subbands. The tunneling 
matrix element is finite only due to the presence of the 
spin-orbit interaction, which allows for spin-up to spin- 
down (and vice versa) conversion. 

The resulting SDW state affects both spin and charge 
properties of the wire. In particular, it suppresses ef- 
fect of (weak) potential impurity, resulting in the inter- 
esting phenomena of negative magnetoresistance in one- 
dimensional setting, as described in Section fVl 

SDW ordering acquires true long-range nature in the 
case of spin chain, where charge fluctuations are absent. 
The staggered moment points along the DM axis and is 
orthogonal to the applied magnetic field. 

Even when the magnetic field and spin-orbital direc- 
tions are not orthogonal, an arrangement when the crit- 
ical Luttinger state survives down to the lowest temper- 
ature (see Section llV Ap . the problem remains interest- 
ing. In this geometry an ESR experiment should reveal 
two separate lines, which represent separate responses of 
right- and left-moving spin fluctuations in the system. 

It is worth pointing out that unusual consequences 
of the interplay of spin-orbit and electron interactions 
are not restricted to one-dimensional systems only. 
We have recently shown^ that Coulomb-coupled two- 
dimensional quantum dots acquire a novel van der Waals- 
like anisotropic interaction between spins of the localized 
electrons. The strength of this Ising interaction is deter- 
mined by the forth power of the Rashba coupling cxr. 

We hope that our work will stimulate experimen- 
tal search and studies of strongly interacting quasi-one- 
dimensional systems with sizable spin-orbital interaction, 
in particular regarding their response to the (both mag- 
nitude and direction) applied magnetic field and/or mag- 
netization. ESR studies of spin chain materials with uni- 
form DM interaction are very desirable as well. 

We would like to thank I. Affleck, L. Balcnts, G. Fi- 
ete, C. Kane, D. Mattis, K. Matveev, E. Mishchenko, J. 
Moore, L. Levitov, J. Orenstein, M. Oshikawa, M. Raikh, 
K. Samokhin, and Y.-S. Wu for useful discussions and 
suggestions at various stages of this work. O.S. thanks 
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APPENDIX A: DERIVATION OF 
INTRA-SUBBAND HAMILTONIAN (p9|) 

Non-interacting (kinetic energy) part of the p-th sub- 
band (v = ±) Hamiltonian reads 

flg.tr. = y Jdx {d x <j> v ? + {dj v f. (Al) 

The intra-subband interaction term is given by the part 
of (|12p which involves only densities from the z^-th sub- 
band 

#intra = \ J dxdx'U{x - x')p v (x)p u {x') , (A2) 

where the density in the u-th subband is expressed with 
the help of (JT^J) as 



Pv 



= rIr v + lIl + cos [7^ 



x e 



-i2k v x 



A2k u x : 



LtRj. (A3) 



Its bosonized form follows 

Pv = ^=d x d> u - sin[V4^^ + 2k v x] , (A4) 

V7T 7ra 

where the first (second) term represents uniform (2k u ) 
parts of density. The interaction term (|A2|) then natu- 
rally splits into a sum of two contributions 



-^intra ~ ^0 + ^2k v > 
r, _ ^(0) 



H' = 



Hoi 



2n 



COS 2 ^] 



(A5) 



dxdX U{x) sin [V 4tt(/v (X + x/2) + 



l 2k v — r,/ \2 

+2k u (X +x/2)} sin[V^r(f> v (X-x/2) + 
+2k v (X - x/2)], 



(A6) 



where x — > x — x' and X — (x + x')/2 are the relative 
and center-of-mass coordinates, respectively. 

Next, following Ref. we fuse the two sines in (|A6j) 
(the result is denoted as S below) using (|26|) and OPE 
identities (|B1IB2|) to obtain 



(A7) 



x exp{i/zV47r[<KA + x/2) - cf>(X - x/2)}}. 

Performing gradient expansion in x, neglecting boundary 
contribution (equivalently, using periodic boundary con- 
ditions so that J dX dx4>{X) = 0), and summing over 



p = ±1 leads to 



H 2k u - 



cos 2 [7^ 



2tt 

X / dX {dx^f 



dx U{x) cos[2fc„ir] 



(A8) 



The integral over relative distance gives backscattcring 
component of the potential U(2k,y). Adding two contri- 
butions we find 

The sum of ([ATI) and (fAT)]) gives us the result ([29]) . 



APPENDIX B: PERTURB ATIVE APPROACH TO 
GENERATE THE COOPER TERM 

Expansion (|73l) relies on the following operator product 
expansion (OPE) 



e ia<p R (z) .. e it34> a (z') 



xe 



- a p(<t, R {z)<l> R {z')) 



xe 



Jla<p L (z)+P<f> L (z')] . 
-al3(4, L (z)<t> L (z')) 



(Bl) 



where z ~ u„t + ix and z = u a r — ix and correlation 
functions of chiral bosons are defined by 



o L (z)<p L {z')) = -J-ln(£_!), 

47T do 



^h(z)M^) = --^m( — ). 



47T 



<l I) 



(B2) 



Both (|B1[) and (|B2|) follow from the harmonic 5*0, see 
(f?4")) . We also employ Baker-Hausdorff formulae 



A p B = p B p A p [A,B] p A p B = p A+B ^[A,B] 



e e 



(B3) 



to convert expressions in terms of dual bosons ip a , 6 a into 
those in terms of chiral bosons 4>R, 0l, 



(B4) 



(Note that for brevity we suppress spin index a on the 
right-hand-side of the above equation.) Their commuta- 
tion relations are given by (|2"3")l and (f2"4")l with v = v' = a. 

Series (|73p are conveniently formulated, using (|B3[) and 
K a — » K, in terms of 

_. p iV2TrKfi(p <r piy^fvSv ■ jMox _ 



A,j, v (z,z) =: e 

= e l f (7?7-^)-*f . e iv^F(/iV / tf+^=)<M*) 



xe 



Indeed, 



H R = I dx ^ vAn, v . 

fJ,,U=± 



(B5) 



(B6) 
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In the second-order (g R ) term, we need to combine 
A^ v (z, z) and A, Jl > v >(z', z'). Using (|BTj) and (|B2]) , we fuse 
field at different points z, z' — u a r' + ix' by set- 

ting z' — > z everywhere where this procedure does not 
cause divergence. This is just an OPE-based gradient 
expansion. In this way we obtain 

A ia ,(z)A lt , v ,(z > ) =ex P [i|(-i-if)-^(/^ + /iV)] 

x exp[i-(fj,fj,'K ~ — + iiv' - vu.')] x (Aio^+n'*') 
2 K 



x : exp 



W27rK(fJ, + fi')<f>L(z) + i 




{v + v')4> L {z)} : 
/iz/ + v/J,')] 



x : exp[iv27rK(^ + p')4> R {z) - 
x exp[-- ln(— - 7 )( MM K +-£T 
x ^[-l^-fo-^'tf + ^- - ^ - „ M ')]. (B7) 

In the subsequent summation over /i, fjf = ± indices 
two combinations, with fi = fjf , v = v' and = — //, v = 
— v', produce highly irrelevant terms (of scaling dimen- 
sion w 4) and are disregarded easily. The choice fi = p! 
and v = —v' produces backscattering correction 



A liV (z)A fl - u (z r ) = exp[iir(K + 1/K - p>v\ 



an 



K-l/K 



(B8) 



Finally, the choice /i 
vant Cooper term, 



-// and v = v 1 , yields the rele- 



ApvA-pv = exp[z7r(//i/ - K - 1/K)] 



«o 



K-l/K 



ip.qo{x-x') iv^/Sn/K 6 a 



(B9) 



We then notice that for A" = K a > 1 (see (gTJ with 
7_f — * 0), the backscattering piece (|B8|) is not singular as 
z' — > z limit is taken and simply disappears in this limit. 
Moreover, it is irrelevant (scaling dimension 2/K a > 2) 
and contains an oscillating phase factor 2qQX which "av- 
erages" it to zero. The other, Cooper contribution (|B9[) 
instead diverges in this limit and thus has to be retained. 

Collecting everything together, we find for the second- 
order correction to the unperturbed 2" (0) = / e" s ° 

Z (2) = f dXdToa&y^/K B a {X,T)], 

8 u a J 

(BIO) 

where (X,T) = ((x + x')/2, (r + t')/2) are the center of 
mass coordinates. Function f{n), with n = (K—l/K)/2, 
is given by the integral over the relative coordinates 
(x, t) -> (x - x', t - t') (it comes from OPE result (|B9])) 

cos[qox] 



/(«) 



(x 2 + i 2 )- 



-2v / 7rcos(7rK) 



r(2 - 2«)r(«- i/2)/r(K), 



/(«-0) = ±?. 

9o 



(Bll) 



Because of its convergence, the integral can be ex- 
tended to infinite limits, and this was done here. Re- 
exponentiating this contribution, we end up with the 
second-order correction SS^ to the free action So 



dXdT cos 




(B12) 



Note that U(2kp)/vp comes from small-K limit of (|B11[) 
while the ratio of spin-orbit to Zeeman energies ap- 
pears from gji/qQ combination. Finally, observe that 
Klein factors combine to produce overall positive sign, 
as (mVl f = mVimVl = -rfvl = -I, in perfect agree- 
ment with our two-subband result in (|37[) . 



APPENDIX C: MOMENTUM SHELL RG 

This appendix is intended to provide self-consistent 
description of the standard momentum-shell RG to the 
problem and to highlight few seemingly tricky technical 
points that arise. 

Similar to the Appendix [B] our starting point here are 
equations (|73I74[) and (fT2l . Fields <p a ,6 a are split into 
slow (index s) and fast (index /) components 



if(k) = , 8(k) = d r(k), 



(CI) 



r=s,f 



i~=s,f 



where k = (q, u>) is a two-momentum. Fast components 
have finite support only in the narrow (two-) momentum 
shell of "width" dA, and are integrated over during the 
RG step. Precise shape of the momentum shell to be in- 
tegrated out will be discussed near the end, most of the 
calculation relies on the fact that the area of that shell 
is small, cx dA. Being quadratic, the free action splits 
into slow So,s and fast Sqj parts as well and, integrating 
out fast modes in every order of the perturbation expan- 
sion, one converts (|73[) into a cumulant expansion for the 
effective action 



1 



SeS = S , s -{S') f --((S' 2 ) f -(S') 2 f ) + 



(C2) 



where the perturbation is S' = J drHpi and (O) / stands 
for the expectation value of O evaluated with fast action 
Sqj. The first order term just produces rescaling of the 
coupling constant 

9R - 9r = ~9R x exp[-i(tf + l/K) f y], (C3) 

where J denotes integration over the fast component 
support area. The factor in front of the integral is 
just the scaling dimension of the spin-orbit operator 
Hr, h(K + l/K) w 1. Rescaling space-time back pro- 
duces additional factor of 2 in the exponent in (|C3|) : 
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\{K + 1/K) -f 2 - \{K + l/K). This, using /' f 
dA/A = d£, leads finally to the standard first-order RG 
equation for the running coupling constant, dgn/dl = 
(2— \(K + 1/K)) (Jr. Note that Klein factors rj^i do not 
affect the scaling in any way. The second-order contribu- 
tion contains Cooper and backscattering terms, 

6S (2) = J J[B(r l2 ) - 1] cosfv^tffov - <p s>2 ) 



-qo(xi - x 2 )\ cosU 



+ 6 sa )]-[B-\ ri 2)-l] 



x cos[v27T^(^ S! i + <p S)2 ) + qo(xi + x 2 )] 
1 2% 

x cosU/ — (0 3 i - e St2 )], 

V iv 



(C4) 



where 1(2) are short-hand notations for r\, 2 — (^1.2, T\ >2 ), 
r*i2 = ?\ — r 2 , and, for example, ip Si i = ip s (xi,r{). Here 

B(r 12 )=eM(K-l/K) f ^^_DA } (C5 ) 



2tt 



k 2 



appears from integrating out fast components of bosonic 
fields. Performing gradient expansion of (|C4j) we find two 
contributions 



SS, 



(2) 



SS 



(9' R ) 2 
8 

-{9'r? 



Br 



d R cos 




B hs / d 2 Rcos[\/8irK(p s (R) 



-2q X] 



(C6) 



(C7) 



where the center-of-mass coordinates are as usual R = 
(X,T) = (fi + r 2 )/2. Their amplitudes are given by 

B c = J d 2 r[B(r)-l]cos[q x], (C8) 

S bs = /" d 2 f [S-^rJ-l], (C9) 

and differ by the absence of cosine factor in the second, 
backscattering-rclatcd amplitude (IC9|) . The final step is 
to expand B(r) as according to its definition (|C5[) the ex- 
pression in the exponential is proportional to small dA 56 . 
One then finds that (|C8p is proportional to the product 
of the integral over the relative coordinate, J d 2 r and 

the integral over fast mode support, J d 2 k, coming from 
(|C5[) . The coordinate integral is performed first to obtain 



n(K-l/K) / diudq 



5(u)[S(q-q )+5(q + q )] 



Be 

ui* + 

(CIO) 

We finally argue that magnetic field, which determines go 
(|72[) . breaks the symmetry between space x and time r 
and this allows us to choose an asymmetric prescription 
for J . Namely, we integrate over all frequencies while 
the g-integration is restricted to the ±|A — A'| interval. 
This gives 



Br 



_ / 2n(K - l/K)/q 2 if g e (A - dA, A) 
otherwise. 



(Cll) 



Note that the result does not contain dl = dA/A 
which serves to emphasize its meaning as a fluctuation- 
degenerated initial value of the Cooper term. Let us now 
see what this approach predicts for the backscattering 
term (IC9[) . Expanding B(r) as in (jClOp we get 



B hs = -2^K-1/K) f dudq 5 -^ 5{q) 
J w 1 

+ l(K-l/K) 2 fdud qw L_ 



2 +q 2 



(C12) 



While the first term is clearly zero, the second, which 
originates from the second order expansion of (|C5p . is 
clearly finite for any shape of the fast modes support. 
This, quadratic in U(2kp)/vF <§C 1 result, is in agreement 
with a non-singular structure of the similar correction 
found during the real-space calculation in (|B8|) . Note 
finally that in the absence of magnetic field (go = 0) 
this scheme predicts similar in structure (but different in 
signs) corrections to the Cooper (|C6[) and backscattering 
(1C7[) terms, in agreement with the result of real-space 
OPE calculations, see Chap. 20 in Ref.[[3(| and papers 

Now we return to the original problem with finite go =/= 
0. Combining (|CTT1l with |C6} we have 



SS, 



(2) 



/ a R k F \ 2 U{2k F ) f 2 - 
\A z kJ 4(7ra ) 2 i " 



cos 




0s(R)], 



4(7ra ) 2 J L V K 

(C13) 

which agrees with (|B12|) and (|37|) in everything but signl 
That sign comes from the Klein factors in g R (note that 
at this stage the difference between c/r and g' R is of higher 
order in d£ and not important) and is a consequence of 
the identity (npn) 2 = —1. This puzzling discrepancy be- 
tween (|37IB12[) and ([CT3)) is worth fi guring out in detail. 

We note that within the functional integral approach, 
which is the framework for the momentum shell RG de- 
scribed here, all information on commutation relations of 
dual fields (p and 9 is contained in (—id T (pd x 9) term in 
the bare action ((74)) . This is nothing but field-theoretic 
version of the canonical pi term in quantum mechanics 59 . 
It identifies ip as a "coordinate" and d x as a "momen- 
tum": [ip(x),d x '9(x')] = iS(x—x'). While fully consistent 
with our basic commutation relation (|68p . this canonical 
bracket does not contain information on the non-trivial 
commutation relation (|2"3")l of chiral right and left bosons. 
Indeed, it is a simple exercise to see that replacing (|23|) 
with commuting chiral bosons [4>r, 4>l] = changes our 
into [<p(x), 0(x')] = i sign(a/ — x)/2 which still satis- 
fies "coordinate-momentum" identification for the pair (p. 
and d x 9. To put things differently, an analog of bosonic 
action ([74")) expressed in terms of "tilded" fields (p and 9 
is identical to the current one in terms of our ip and 9 
fields satisfying equations (|2"31 12"4"1) and ([55]) . This simply 
means: bosonic action Sq, (|T4[) . does not enforce anticom- 
mutation relations between right R and left L moving 
fermions (|22p . This shortcoming of bosonic functional 
integral is well known, see for example Appendix C in 
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Ref. |3jJ, and several "fixes" were proposed in the litera- 
ture. Our approach consists in enlarging the role of Klein 
factors: instead of (l22l) we bosonize fermions here as 



R K = 



KRs 



, L s = ^^ e - lV4 ^^.(C14) 

V27ra \ / 2Tra 

The Klein factors kr/l,s now carry double index: chi- 
rality (R or L) and spin (s). Even though the chi- 
ral bosons (f>R s and 4>l s now commute, the anticommu- 
tation of R s and L s is enforced by the Klein factors: 
{k\, s , Kyy} = 25\ : yS S:S >, with X — R/L. 

As a result of the proposed modification the spin-orbit 
ter m (172")) has to be modified. It is convenient, following 
Ref. 60, to introduce 



(C15) 



This product satisfies f 2 = 1, from where it follows 
that its eigenvalues are T = ±1. One also checks that 
[f, Kr^Kri] = and kr^kriT = — kl^kli- These prop- 
erties allow us to represent Hr as 



Hr = Gr J dx{ ^ * T ' cos[V27r^ 



q x\ smU/ —0\ 
V A 



- (1 - F) sin[v^gy + gox] cos[ifeg]}, (C16) 
2 V A 



2a R k F , 

<^fl. = (ikrtKri, 

7ra 



(C17) 



Finally one observes that f commutes with the Hamil- 
tonian of the problem, which implies that its eigenvalues 
represent integrals of motion. That is, the choice T = +1 
or r = — 1 is the gauge choice, and one can replace op- 
erator r in (|C16j) by its eigenvalue T. Comparing (|C16j) 
with our original (|72[) we observe that the latter corre- 
sponds to r = +1 gauge: with T = +1, equation (|C16|) 
transforms into (fT2")) by replacing (tp,8) — > (<p ai a ). Re- 
peating steps that led us to (|C13[) we arrive at 



SS, 



5(2) 



4 V qa ) 



l/K) I d 2 R cos 




(C18) 



The only, but key, difference with (|C13j) is that now 
(iK R ^K Rl ) 2 = -(kr^k ri ) 2 = K 2 R1 K 2 Rl = +1, which im- 
plies that the amplitude of the Cooper term is positive in 
(|C18[) , in final agreement with our previous and indepen- 
dently derived results in (|37IB12[) . This amusing exercise 
illustrates the importance of Klein factors, and, more 
generally, of preserving correct (anti)commutation rela- 
tions when implementing convenient but tricky bosoniza- 
tion formalism. We conclude by noting that results of the 
other gauge choice, T = — 1 in (|C16[) . while equivalent in 
principle to the one made above, are most conveniently 
understood as following from the global shift of bosonic 
fields: ip -► ip+y/n/(8K) and 6 -> O+y^W/S. This shift 
must be made in all bosonized expressions: it changes 
overall sign in (|C18[) but this is "compensated" by the 
effect of the global shift described here. As a result, the 
physical meaning of the Cooper instability as that of the 
SDW X instability remains intact. 



APPENDIX D: PERTURB ATIVE APPROACH TO 
GENERATE THE IMPURITY TERM 

The purpose of this section is to show that inter- 
subband contribution to impurity potential, second term 
in Eq. [86l can be obtained as a result of interference be- 
tween local impurity backscattering (first term in (|86p ) 
and bulk spin-orbit (|70p . Being interested in the interfer- 
ence between two single-particle terms, we perform the 
calculation directly in fermion fields and specialize to the 
limit A z 3> E s _ (jp — > 0) for simplicity. 

The lowest order in the perturbative expansion involv- 
ing these two terms is obtained by the following correc- 
tion to partition function (compare with l|73p) 



3Zi mp — — 



\ J e- S " J dTdT'H R (T)V B (T') 

-» J e- So J dr SV B , (Dl) 

where the second line identifies correction to the impurity 
backscattering term we are after. Here Sq is the action of 
two Zeeman-split {T>4} subbands, and V B is the intra- 
subband term due to the impurity, located at i = 0, 
given by 



V B = V 



(//;/- + L\R t 



a;=0 



(D2) 



The spin-orbit term Hr in the presence of magnetic field 
(so that fcf — ki = 5k f) reads 



Hr = a R k F I dx V ( R\R- K e- ls5kFX - ^ T - is5kFX 



(D3) 

where we neglected terms oc Skp/kp = A z /Ep <C 1. 

We calculate ([Dip by fusing fermions with like 
spin index, that is by making the replacement 
^r/Ls(x,t)^(x' ,t') -> G R / L (x- x',t-t') where pos- 
sible. Here Gr/l stands for fermions Green's functions. 
To lowest (zeros) order in the interaction these are given 
by the free fermion Green's functions 41 



Gr(x,t) = (* Rs (x,t)* Rs ,(0)) = 
G l (x,t) = (* Ls (x,r)*L,(0)) = 



2it(vt — ix) 

$s,s' 

2tt(vt + ix) 



(D4) 



In this way we find 



dr SV B = aRkFV ° 



J dr^lt-L^P, 



where 



I 

Pi = J dxdt( 
P 2 = dxdt( 



4tt 

-{r]l ] -l]r ] )p 2 



-—iSkp-x 



Z 

-iSkpx 



(D5) 

(D6) 
(D7) 
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and z = vpt + ix, z = 
are easily calculated 



VFt — ix as usual. These integrals 



Pi = -Pi = 



An 



SkpVp 



(D8) 



Expressing Skp = A z /vf we obtain the correction 



5V M = 



(i?jL i -i?{i T +L t i i? T -L t T i? i ). (D9) 



We observe that SV B is odd under spatial inversion V 
(with respect to impurity location) when x — > —x and 
right- and left-movers get interchanged, R s <-> L s . Of 
course, it must be odd under V as it is obtained from 
fusing even (|D2[) and odd (|D3|) in V terms. The oddness 



of (|D9|) is the reason for the relative minus signs in this 
equation. 

Bosonization of (|D9|) . following Sec. IIII Al results in 



SV 1 



2a R k F V irjirji 



ira 



^[V^<Pp]cm[V2TrO<r], (D10) 



which confirms our previous result (|85l86p . The gener- 
ated impurity potential describes intra-subband impurity 
backscattering. 

The absence of the potentially dangerous term with 
sin[\/27r0<T] in place of cosfv^^tr] in (|D10[) is now clear: 
such a term would require no minus signs in (|D9[) which 
is forbidden by the oddness of 8V B under inversion V . 
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